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Abstract 

We have investigated the correlation functions of interacting bosons at the generic superfluid- 
insulator transition, a prototypical quantum phase transition, in two dimensions in the spherical 
limit. Unexpectedly the spatial correlation functions show non-power-law behavior consisting of 
two parts: short-range correlation due to the particle-hole pair excitations and long-range off- 
diagonal order due to the single-particle condensation. The temporal correlation functions, on the 
other hand, show power-law behavior. 
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Quantum phase transitions of interacting bosons 2] have continuously drawn significant 
attention. They have been realized in many different systems, such as Josephson-j unction 
arrays^, thin film superconductors 4], 4 He films)^] and ultracold atoms in optical lattices^, 
Q|. Typically the transition occurs between a superfluid and an insulator, and often serves 
as a prototype of the quantum phase transition[8( because of its convenience in describing 
the transition in terms of the establishment of a macroscopic phase coherence. 

Quantum phase transitions are frequently discussed in equivalent classical models which 
include fluctuations in the temporal direction. For the special case with the dynamical crit- 
ical exponent z = 1, the equivalence of the spatial and temporal directions directly confirms 
the validity of this mapping. In this scheme the properties of quantum phase transitions 
in d dimensions are readily derived from the classical models in d + 1 dimensions. In gen- 
eral, however, asymmetry between the spatial and temporal directions in quantum phase 
transitions results in z ^ 1. Even for this case, a simple extension of the above scheme 
into d + z dimensions is frequently assumed. Many scaling ansatz, such as hyperscaling 
relations [sj], are constructed on the basis of this assumption. Similarly, correlation functions 
at the transition are assumed to have power-law behavior effectively in d + z dimensions. 
Some numerical works have used this assumed property to extract the value of z from cor- 
relation functions h], U|- However, it is quite plausible that the asymmetry of fluctuations 
in the spatial and the temporal directions, which brings non-unity of z, would modify the 
correlation functions. 

In this work, we investigate the scaling properties of the quantum rotor model, which 
is equivalent to the boson Hubbard model in large density limit, at the generic superfluid- 
insulator transition in two dimensions in the spherical limit. We find that the scaling behav- 
ior clearly supports hyperscaling in d+z dimensions with z = 2, as predicted^]. However, the 
spatial correlation functions at the transition show, instead of power-law behavior, properties 
consisting of two parts: one is a short-range contribution from the particle-hole pair excita- 
tions and the other is the off-diagonal long-range order due to single-particle condensation. 
This implies that the macroscopic phase coherence is achieved through the single-particle 
condensation rather than the establishment of long-range correlation of fluctuations in a 
spatial direction, and we have the background particle-hole excitations as a normal fluid, 
instead of the single-particle excitations, even in the case with a charge-offset. The temporal 
correlation functions, on the other hand, show the power-law behavior as expected. 



2 



The universal features of quantum phase transitions in interacting Bose systems are di- 
vided into two universality classes 3] depending on whether the transition is driven either by 
the phase fluctuations or by the density fluctuations of bosons. At a commensurate density 
the particle-hole symmetry is sustained across the transition, and the transition is driven by 
the phase fluctuations or, equivalently, by the particle-hole pair excitations. This transition 
belongs to the (d+l)-dimensional classical XY spin model with z = 1. In the presence of a 
charge-offset, however, either single particle or hole excitations are favored: the transition 
is marked by the disappearance of the energy gap for single particle or hole excitations. 
Therefore, the universality class of this transition, known as the generic superfluid-insulator 
transition, is characterized by the single-particle nature, carrying z = 2 and the mean-field 
correlation length critical exponent v — 1/2. It is interesting to investigate whether corre- 
lation functions could reveal directly the different nature of fluctuations in two universality 
classes. 

The essential properties of strongly correlated interacting bosons can be captured by a 
boson Hubbard model 

#=f]p»?-pX>-*E ( 6 & + b H c 1 ) 

i i <ij> 

where bj(bj) is the boson annihilation(creation) operator at j-th site, and rij is the number 
operator. U and t stand for the strength of on-site repulsion and nearest neighbor hopping 
respectively, and u is the chemical potential. It is convenient to put \ij\J — > no + u/U with 
an integer no to make —1/2 < u/U < 1/2. Here no is an integer representing the background 
number of bosons per site. We study the phase transition of this model on two-dimensional 
square lattices. 

The phase transition of this model is characterized by establishment of the phase co- 
herence of the order parameter. In the limit n ^> 1, the Hamiltonian is reduced to the 
quantum rotor model 

E = \ E^-^E n *- 2J E cos ^-^-)> (2) 

where 9j is the phase angle of the order parameter, nj = —id/dOj, and J = not. 

Now we use the spherical approximation to investigate the quantum phase transition of 
the above model. The spherical approximation jl2j has been widely used in studies on phase 
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transitions, which is very powerful in handling models beyond the mean-field level. This 
method treats phase transitions in a scheme which is exact in the limit that the number 
of order parameter components goes to infinity. It has been used to study quantum phase 



transitions in Bose systems 13| |. Previously this method was used lj] to study the superfluid- 
insulator quantum transition at a commensurate density in two dimensions, yielding v = 1 
and z = 1, as expected. 

Through a path integral mapping, we can construct the corresponding classical action 

sbM = - / ^[tfJ^toCar-AO 2 ^ (3) 

where K = 1/(2U), (3 denotes the inverse temperature, denotes the nearest neighbor 
hopping matrix elements, and ipj — e ^ l9] ■ The trace in the partition function should be 
taken with the constraint = 1 for every j. But we introduce a Lagrange multiplier a 
to replace the constraint by a less strict self-consistent condition 1 = (^j(r)^*(r)) , where 
(...)o represents the average over So, which is called the spherical approximation. 

In this work we focus on the case with a charge-offset (n = fi/U ^ 0). Scaling theories 
predicted that the dynamical critical exponent z = 2|2]. To confirm this through the hyper- 
scaling hypothesis in d + z dimensions, we investigate the finite-size scaling behavior of the 
superfluid stiffness 

p s = L^ d +^X p (L^(K - K C )J/L Z ), (4) 

where X p is a scaling function. By diagonalizing the equation ^(^i — Jij)^ = to 
have (f>l = (1/ \fl?')e % ^ Ri and e^* = a — 2 J(cos q x + cos q y ), where Rj is the position vector of 
j-th site and q = (2ir/L)(n x , n y ) is the wave vector with integer n x and n y , we can easily 
derive the formula 

P' = J2 E {^WT^ C °* q * - ^r*) [ooth ^(<sf*/K- fi) + coth ^{sJJ/K + /, 
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J 2 psin 2 q x r 1 + 1 O 



L sinh 2 f(y^7^-/x) sinh 2 f(^7^ + /i)- 
Here we set the lattice constant to be 1, and take the energy unit J — 1. 

Figure n shows the finite-size scaling behavior of p s for n = 0.1. With z = 2, curves 
for different sizes L = 60, 80, 100, 200 cross at a point K c = 0.09513(2). Here we fix the 
aspect ratio to have (3/L z = 0.02, and obtain the correlation critical exponent v = 0.5, as 




0.09508 0.0951 0.09512 0.09514 0.09516 
K 



FIG. 1: The finite-size scaling behavior of the superfluid stiffness, p s for n = 0.1. High quality 
crossing behavior at K c = 0.09513(2) with z = 2 for L = 60, 80, 100, 200 is obtained without 
suffering from any statistical errors. We set (3/L z = 0.02. Inset: The data collapse onto a single 
curve as a function of the scaling variable L l / U {K — K c ) with v = 0.5. 

shown in the inset. The high quality crossing behavior implies that p s D- d+z ~ 2 ' is indeed a 
universal quantity in the vicinity of the critical point, strongly supporting the hyperscaling 
hypothesis in d + z dimensions. 

Now we turn our attention to the correlation functions. The correlation function between 
Ri and Rj with distance r in time is given by 
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G(|A-^|,r) = I? X) 



•[" - al + 



-] 



(6) 



2\IK<3 L i_ e -/3(V^+^ 1 - e -«V^-f) 

for < r < 1. It was assumed fl Q that at the critical poiet the co rrela t 10 n tactions 
show long-range power-law behavior in d + z dimensions so that G(\Ri — Rj\, r) ~ 1/(|-Rj — 
Ri\ 2 + r 2 / z W 2 , where w = <i + z — 2 + r? (r? is the anomalous critical exponent). Based on 

nnn 

this assumption, estimation of the value of z has been made[10, by comparing the 

asymptotic behavior of the spatial and temporal correlation functions. 

Figure El shows the spatial correlation functions along a spatial direction at the criti- 
cal point. Surprisingly they do not show the power-law behavior; instead, the correlation 
functions have forms 



G(x, 0) = A(x- l e~ x ' x + (L — x)~^e 



c 
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FIG. 2: Log- log plot of the correlation functions along a spatial axis at the critical point K c = 
0.095132 for n = 0.1. They do not show the power-law behavior. The curves are fitted perfectly 
onto the form G(x, 0) = A{x~ x e~ x l x + (L - x^e'^'^^) + C. 

for 1 <C x <C L, where A is a finite number, insensitive to the size of the system, C is a 
constant depending on the size, and A is a fitting parameter. For example, for the curve in 
Fig. H with L = 200, A = 0.26, A = 6.11, and C = 1.17 x 10" 4 . Similar features found in 
Monte Carlo calculations were discussed as a result of finite-size effect j^5 1 . Here we use the 
same value of the aspect ratio f3/L z = 0.02, previously used in Fig. ^ For significantly wide 
range of /3/L z , the same behavior is found. Very small aspect ratios, however, change the 
behavior, possibly due to finite-temperature effects, but still no power-law behavior occurs. 

We can easily identify the meaning of the constant C. The spatial correlation function 
is actually flat for x ~ L/2. In other words, C = G(L/2,0) = \ (ip)\ 2 , where 1(^)1 is the off- 
diagonal long-range order parameter. Note that the constant C is absent for the case n = 
in which the correlation functions show pure power-law behavior G(x,0) ~ x _1 jlJ . Ilc| . 
Figure H3 shows the scaling behavior of We obtain the mean-field critical exponents 

v = 0.5 and /3/v — 1, where f3 is the critical exponent characterizing the scaling properties 
of the order parameter. 

The long-range off-diagonal order is coming from the single-particle condensation. We 
can find this clearly by observing the momentum distribution, n q , in Figure EJ Here we 
plot the momentum distribution as a function of q = (2ir/ L)(m, m). The density of the 
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FIG. 3: The finite-size scaling behavior of the off-diagonal long-range order parameter for L = 
100, 120, 200, providing the mean-field critical exponents v = 0.5 and (3/v = 1. 




FIG. 4: Momentum distribution of the quantum rotor model near the superfluid-insulator transi- 
tion. The density of condensate at q = rises sharply as the superfluid onset transition is tuned. 

single-particle condensate at q = state rises as the superfluid onset transition is tuned. 

The momentum distribution consists of the single-particle condensate at q = and back- 
ground normal fluid, bringing the macroscopic phase coherence and contributing to short- 
range correlations, respectively. It is quite interesting to note that the normal fluid of 
interacting bosons is not the single-particle excitations but the particle-hole pair excita- 
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FIG. 5: The range of the correlation due to the particle- hole excitations as a function of the 
charge-offset n. It diverges like A ~ 0.61/n as n —* 0. 

tions which contribute to the asymptotic correlation G(x, 0) ~ x~ v e~ x ' x with y = 1 rather 
than y = 2. The correlation range of the particle-hole pair fluctuations at the transition 
diverges like A ~ 0.61/n as n — > 0, as shown in Figure El ultimately yielding the power-law 
correlations at n = 0. 

The temporal correlation functions, however, show the power-law behavior G(0, r) ~ 
l/ T ( z +v)/z w ith [z + rj)/z = 1.0, implying i] rs 0.0. This exponent and other exponents 
obtained above confirm the hyperscaling relation 2(3 = v(d + z — 2 + if). But it is not 
possible to determine the value of z directly from the correlation functions if they do not 
show the power-law behavior. 

In this work, we have studied the generic superfluid-insulator transition in interacting 
Bose systems via the spherical approximation. Scaling properties of physical quantities, 
such as the superfluid stiffness and the off-diagonal order, confirm the hyperscaling relations 
in d + z dimensions, providing the critical exponent z = 2, v — 1/2, and (3 = 1/2. These 
values of the critical exponents, supporting the mean-field nature of the transition, are con- 
sistent with theoretical predictions and other calculations. However, the spatial correlation 
functions at the transition do not show the power-law behavior. Instead they indicate that 
the correlations consist of two parts: short-range correlations due to the normal fluid of 
the particle-hole pair excitations and macroscopic long-range off-diagonal order due to the 
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FIG. 6: The temporal correlation functions at K = 0.095132 for n = 0.1. They show the power-law 
behavior G(0, r) ~ i/ T ( z + r i)/ z with (z + rf)/z = 1.0. It implies that r/ ~ 0.0, but the determination 
of z is not possible. 

single-particle condensations. The temporal correlation function, on the other hand, show 
the power-law behavior even though direct determination of z is not possible. 

This work was supported by grant No. R05-2004-000-1 1004-0 from Korean Ministry of 
Science & Technology. 
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